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Abstract 



We suggest a new, renormalization group (RG) based, nonperturbative 
method for treating the intermittency problem of fully developed turbulence 
which also includes the effects of a finite boundary of the turbulent flow. The 
key idea is not to try to construct an elimination procedure based on some 
assumed statistical distribution, but to make an ansatz for possible RG trans- 
formations and to pose constraints upon those, which guarantee the invariance 
of the nonlinear term in the Navier-Stokes equation, the invariance of the en- 
ergy dissipation, and other basic properties of the velocity field. The role of 
length scales is taken to be inverse to that in the theory of critical phenomena; 
thus possible intermittency corrections are connected with the outer length 
scale. Depending on the specific type of flow, we find different sets of admis- 
sible transformations with distinct scaling behavior: for the often considered 
infinite, isotropic, and homogeneous system K41 scaling is enforced, but for 
the more realistic plane Couette geometry no restrictions on intermittency 
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exponents were obtained so far. 
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I. INTRODUCTION 



A fundamental, yet unsolved problem of fully developed turbulence is the scaling be- 
haviour of the structure functions in the inertial range for Reynolds numbers approaching 
infinity. 

Although many if not most experiments seem to favour the existence of intermittency 
corrections [Q^J, some other experiments seem to confirm the classical "K41" exponents 
|J|7]]. Thus, the possibility remains that measured intermittency corrections are due to still- 
to-small Reynolds numbers, invalidity of Taylor's hypothesis, lack of statistics or other 
experimental deficiencies. 

On the theoretical side, there hasn't been much progress in this regard either since 
K41 theory [p|-|T^1 and mean-field theories such as Ref. |13| and Kraichnan's LHDIA |14 



Numerical simulations |T5|-[i7| still do not achieve high enough Reynolds numbers, and ap- 



proximations allowing higher Reynolds numbers such as REWA ||T8|-pO|] reduce the degrees of 
freedom without control of the consequences. Phenomenological models pT| |33| are not only 



limited to physically "explaining" intermittency corrections without being able to connect 
intermittency with the Navier-Stokes dynamics, but cannot even demonstrate its existence 
for Navier-Stokes turbulence since statistical models generically show intermittency anyhow 



The basic motivation to analyze turbulence with RG methods comes from a couple of 
similarities with critical phenomena: the importance of extremely many spatial degrees of 
freedom, strong fluctuations, and nonlinear interaction, to name a few. These ingredients 
are essential for spatially extended self-similar structures and scaling laws with nontrivial 
exponents. We now briefly discuss some of the many efforts to derive inertial range scaling 
exponents with the aid of RG and/or field theory; current overviews can be found in Ref. 



The first, by now "classical" approach of Ref. p9[-[4"l"H formally adopted the formalism of 



Wilson style respectively field-theoretic RG theory for critical phenomena to hydrodynamics. 
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Microscopic degrees of freedom are eliminated by averaging over an artificial, stochastic, and 
microscopic volume forcing. The intermittency exponents depend on the forcing's spectral 
strength, and these authors explicitly state that such a theory is not adequate for fully 
developed turbulence. Nevertheless there are several extensions of these theories such as 
Ref. pE2| , pE5|| where the incorporation of traditional closure methods enables the calculation 
of dimensionless numbers like Kolmogorov's constant. Other work employs a driving force 
restricted to small wave numbers without being able to calculate an effective, RG invariant 



forcing ||44]| . Yet another idea is to examine freely decaying turbulence with random initial 
conditions |j45| . Last but not least, formal analogies to other RG theories such as renormalized 
quantum field theory are possible |37| . 



Quite early in the game, another mapping to critical phenomena was suggested which 
seems to match the common phenomenological point of view of turbulence much better 
p6^p0[. It is "inverse" in the sense that macroscopic degrees of freedom should be elimi- 
nated, and that, consequently, intermittency corrections are related to the outer length scale. 
Another, within this picture obvious, but actually rarely suggested idea is to take the energy 
dissipation rate to be the RG invariant quantity instead of the free energy which does not 
seem to have any special role in turbulence P5|JoT| . We are aware of two actual attempts 
for an inverse RG theory, one with Wilson-style elimination of degrees of freedom |48j and 



one field-theoretic one |g2|. Although being inverse, both adopt major basic properties of 
thermodynamic systems which do not seem to be adequate for turbulence. 

All methods mentioned so far use perturbation theory with the strength of the nonlinear 
interaction term of the Navier-Stokes equation being treated as a formal expansion param- 
eter which is problematic without a suitable resummation procedure. Furthermore there are 
unphysical infrared divergencies; a problem which can be solved by using a formal perturba- 



tion approach with velocity differences instead of Eulerian velocities ||53|| . This perturbation 



theory is still being worked on |54|-|56f, but so far not with RG methods. 

Without judging all these efforts, it is fair to say that the calculation of inertial range 
scaling exponents from the Navier-Stokes equation has not been achieved yet. 



There is a very important aspect of turbulence which is mostly neglected in turbulent the- 
ories: the boundaries of the system. Numerical simulations usually employ periodic boundary 
conditions and analytical theories try to take care of the finite system size with a formal 
"outer length" parameter. Energy input, which is required to balance energy dissipation, is 
realized with a volume forcing which might be unphysical. 

In this work now, we introduce a new, analytical and inverse RG method which is based 
on the Navier-Stokes equation and tries to take into account the boundaries. The paper 
is organized as follows: in Sect. || we derive the equation of motion for certain velocity 
differences, which will be our order parameter. In Sect. |T| we explain the physical ideas of 
the new RG approach before actually calculating RG transformations in Sect. |V|- |VI]] . We 



end with a summary and discussion in Sect. VTLl 



II. EQUATIONS OF MOTION 

Our starting point is the Navier-Stokes equation 

d t u + u- V x u + V x p = v A x u (1) 

for the Eulerian velocities u(x, t) and the kinematic pressure p(x, t) of incompressible motion 
of a fluid with kinematic viscosity v in a three-dimensional flow volume Q. 

To be more definite, we consider plane shear flow between two infinite plates (Fig. [[]). 
We define the plates to be located at x 2 = moving with velocity Ul in positive ^-direction 
respectively resting at x 2 = L. Reynolds' number is defined as Re = U^Ljv. 



A. Velocity increments 

The choice of the order parameter is crucial for any RG theory. It is therefore desirable 
to find a closed equation of motion for velocity increments instead of considering the Navier- 
Stokes equation for the Eulerian velocities, since such variables are commonly believed to be 
capable of exhibiting universal behavior. This can be achieved by means of a transformation 

5 



introduced by Monin |57|J58[] . Since it does not seem to be widely known in spite of its 



reintroduction in ||53|| , we present it here relatively detailed. The idea is to measure the 
Eulerian velocity field u in a frame of reference defined by an arbitrary but fixed marker 
particle passively advected with the flow (Fig. g). The Lagrangian trajectory s(t|x ,to) for 
such a marker which is at position x = x for t = t is given by 

s(t|x ,t ) = x + / dru(s(r|x ,t ),T)- (2) 
J t 

Next, we define a "velocity increment" 

w(r, t|xo, t ) = u(s(i|xo, t ) + r, t) - u(s(t|x , t ),t). (3) 

The field w(r, t) describes a velocity difference between two space points separated by the 
increment r, for all times t, but has Eulerian as well as Lagrangian characteristics. Monin's 
transformation from u to w thus takes care of the very special property of flow fields being 
coupled to the actual movement of fluid elements. 
Plugging (|3[) into ([TJ) gives a closed equation for w, 

d t w + w- V r w + V x p = v A r w — s. (4) 

There is no explicitly u dependent term due to Galilean invariance. The term — s describes 
a fictitious force (per mass) because the frame of reference, i.e., the marker, moves nonuni- 
formly. Since it is independent of r, it may be eliminated by taking the difference of (|]) for 
arbitrary r as well as for r = 0. Formally, we introduce an operator S, 

S[q](r,f)=q(r,f)-q(0,t), (5) 

and apply it to (|): 

d t w + w- V r w + S[V x p] = uS[A r w). (6) 

The boundary conditions cannot be expressed in terms of w and r alone since the bound- 
aries brake Galilean invariance. The velocity increments w are defined on a flow volume 
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moving relative to the marker, 



fi(s(t|x , t )) = {r = x - s(t|x , t ) |x G fi} 



with the corresponding surface <9f2(s(t|x , to)). So far our introduction of Monin's transfor- 
mation. 

We now specify how statistical averages are defined. First, a local x-centered space av- 
erage over an a-cube is taken for every space point x. Second, an infinite time average is 
performed. Finally, the limit of the averaging volume approaching zero is taken. We as- 
sume that the properties of the nonlinear fluid dynamics guarantee that all mean quantities 
considered in this work exist and are independent of x 1 , x 3 and t . More formally, we define: 



The second equality follows from det V xo s = 1 which holds because V s u = and Lagrangian 
trajectories are unique. 

From now on we will often use the shorthand notation s(t) and w(r, t) for s(t|x , t ) and 
w(r, t|x , t ), respectively. 



The statistics of the velocity increment field w are nonuniversal in general. Not even 
the mean value (w) = XJ(x2 + r 2 ) — U(x 2 ) will vanish because of the anisotropy inherent 
to the turbulent profile Ui(x 2 ) = (ui) (x 2 ) = U(x2)S ij3 which reflects carries the boundary 
anisotropy. It seems therefore reasonable to subtract the turbulent profile from both Eulerian 
velocities appearing in the definition of w. We call the resulting quantity 




— /(r,r|x ,t ). 



(7) 



B. r-velocity fluctuations 



v(r,t) 



w(r,t)-[U(s(t)+r)-U(s(t))] 



(8) 



[u(s(t) + r, t) - U(s(t) + r)] - [u(s(t), t) - U(s(t))] 
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the "r-velocity fluctuation" . By definition we have (v) = and, of course, v is divergence free. 
The boundary condition for v(r, t) is given by the marker's negative r-velocity fluctuations, 

v ( r ^)L^(s( t) ) = -[ u ( s W^)- u ( s W)]' ( 9 ) 

which, remarkably, is independent of r. The quantity v is expected to have universal prop- 
erties and seems to be the most natural choice for the "order parameter" of turbulence. It 
vanishes for laminar flows. 

The pressure p is eliminated as usual by formally solving a Poisson equation (we employ 
standard summation conventions): 

p(x,f)=p f (x,i)- / d 3 £Gn(x,0%%M£,*K(^)), 
Jn 

where the Green function Gn is determined by the geometry, i.e., d Xi d Xi Gn(x, = 5(x — 
£), Gq(x, \ xedQ = 0. Here, the pressure is assumed to be constant at the closed parts 
of the boundary, otherwise there would be an extra surface contribution to p. While this 
is physically reasonable, it should be remarked that the admissibility of such boundary 
conditions is mathematically nontrivial, since the compatibility condition v A x u = V x p has 
to be fulfilled also at the boundary. For the geometry of plane shear flow one can easily 
verify the integral representation 

1 f°° 

G n (x,£) = -— J dkJ (ky/(x 1 -^ + (x 3 -Cs) 2 )- 

coshQQa - £ 2 | — L) — cosh(A:(a; 2 + £ 2 - L)) 

sinh(A;L) 1 j 

for the Green function. The expression for the pressure may be combined with the nonlinear 
term of the Navier-Stokes equation, 

u- V x u + V x p = V[u- V x u] + V x p f , 

where the operator V projects to divergence free functions: 

P i [q](x,t)=g i (x,t)- / <l\il Cvix.^ih qM.I). 

Jn 
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Finally, we obtain an equation of motion for v by inserting @ into (|6|) and employing 
V r w = V r v + V r U(s + r): 

dtv + SV[\-V r v] = uS[A r v]+F[v]. (11) 

The projector V is defined analogously to V: 

A[q](r,*) = ft (r,t)- / d 3 p d ri G n {r + s(t), p + s(t)) 0„.fc(p, t). (12) 

Jn(s(t)) 

The terms containing the turbulent profile U can be interpreted as driving force (per mass) 
for the r-velocity fluctuation: 

JF[v] = -5D[(U(s + r) - U(s)) • V r v(r, t) + v(r, t) ■ V r U(s + r)] 

(13) 

+ i/5[A r U(s + r)] -5V r pf-u(s,t)-5[V r (U(s + r) -U(s))]. 

The forcing term depending explicitly on u is equal to d t [U(s(t) + r) — U(s(£))]. Such explicit 
u dependence is no surprise since the Galilean invariance, which Monin's transformation 
relies on, is destroyed by the boundary. Thus, ([□]) alone is not equivalent to (|l|) because the 
trajectory s(t) is needed which itself is determined by (|l|). 

We consider ( ITT]) as an equation of motion for the order parameter v with given profile 
U and marker trajectories s(t). In this sense it is analogous to the various time- dependent 
Ginzburg-Landau equations used for RG treatment of dynamical critical phenomena, but 
as opposed to the latter the forcing here is deterministic. No statistical distribution has to 
be introduced as it is usually done for RG theories of turbulence. 

III. A NEW RG METHOD 

In this section we outline our physical understanding how a RG theory could be applied 
to fully developed turbulence. The point of view presented in Sect. |III A| will guide us in 
Sect, m B| where the ideas for a new RG method are developed. 
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A. Inverse RG 



After choosing the r- velocity fluctuation v as the order parameter in the previous section, 
a mapping of quantities of turbulence to those of critical phenomena is given. It is based on 
the idea of an inverse RG transformation pursued in earlier work [f46H 50| . 



The large spatial scales, which are governed by a forcing depending on the flow geometry, 
are considered as nonuniversal in the RG sense, and thus correspond to the microscopic 
interaction scales (or lattice constant) in critical phenomena. The small spatial scales where 
viscosity dominates are seen as universal in the RG sense and analogous to the scales in 
critical phenomena set by the macroscopic system size. In this picture one expects anomalous 
scaling for r above i in analogy to r below the correlation length in the critical case and 
regular behaviour for r — > as for the wavenumber approaching zero in the critical case. 
Thus, the (longitudinal) structure functions 

^)(r) = ((v(r,t)-") m ) (14) 

should scale with non-K41 intermittency exponents in the inertial range, and be analytic in 
the viscous range. We do not study correlations in time and implicitly allege that turbulence 
has a "static limit", i.e., a description of equal-time correlations independent of correlations 
in time. 

At this point one might ask which turbulent quantity corresponds to the correlation 
length. If we take the exchange of the role of large and small scales serious, this ought to be 
a viscous length which we will call I in the following. Since there is strong indication that 
the energy spectrum decays in the viscous range according to the universal law 

E(k) oc (|v(k, t)\ 2 ) oc exp(— £ kj (2ir)) k 13 for 2n/k < t 



with some exponent (3 |[32|| , we may define a dissipation length t{y) through t = 
— limfc^oo(hi E(k))/ k. This is completely analogous to the definition of the correlation length 
through the decay of the order parameter correlation function for spatial distances far beyond 
the correlation length. 
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Because there is no evidence that £ vanishes for a finite value of v, v — > is the "tur- 
bulent limit" corresponding to the critical limit of the temperature approaching the critical 
temperature. Using the language of phase transitions, a turbulent system is always in the 
"disordered phase" because of (v) = 0. Since an ordered phase would appear for negative, 
but physically forbidden v only, we do not expect a phenomenon analogous to spontaneous 
symmetry breaking. As the critical limit is defined for fixed microscopic parameters, we de- 
fine the turbulent limit for fixed macroscopic quantities, i.e., fixed geometry and outer scale 
L as well as fixed outer velocity Ul- This is the precise definition of what we understand as 
u Re — > oo" . Thus, v is supposed to be the only critical parameter as opposed to temperature 
and external magnetic field in critical phenomena; U is analogous to the coupling parameter. 

Recognizing this analogy to critical phenomena, the first step of a RG transformation 
for turbulence should eliminate degrees of freedom on the macroscopic scale L, resulting in 
a new field U which is defined in a smaller system with size L = L/b where b > 1 is the 
usual RG rescaling parameter. In the second step the system is rescaled with b to its original 
size V = bL = L and the equation of motion should take on the same functional form with 
renormalized turbulent velocity fluctuation v', viscosity u', and turbulent profile U', which is 
analogous to the renormalization of the order parameter, the temperature, and the coupling 
parameters in critical phenomena. U should approach a fixed point under renormalization. 
The spatial rescaling implies an increase by a factor of b of all physical microscopic length 
scales, in particular £' = b£ for the dissipation length. 

Coming back to the central question of intermittency corrections, we give the functional 
form structure functions are expected to show for v 

f& m (^r) m/3 (£)* 
D (m \r,u,L)=! 

\D { ™\v,L)r m 

Here, the energy dissipation rate e* is supposed to be finite for v — > 0; this will be defined and 
discussed later. The constants b m are also supposed to be finite for v — > 0. Thus, intermit- 
tency corrections are connected to the outer length scale L, whereas in critical phenomena 

11 



0: 



for £ < r < L 



for r < £ < L 



(15) 



the critical exponents' deviations from their Gaussian values show up as a dependence of 
correlation function amplitudes on the microscopic scale. 

The form is also compatible with the concavity of the curve 5( m as a function of m, 
as a simple application of Schwarz' inequality applied to (v m ) shows, leading to 2 5( mi+m2 > 
<5C2mi + ^C2m 2 if L defines the correction scale, so r/L < 1 P5| , |55| . But this argument does 
not rule out more complicated scenarios such as the simultaneous occurrence of L and i as 
correction scales. Also, we have assumed that structure functions are isotropic in the inertial 
and viscous ranges despite of an anisotropic geometry and that they are independent of X2 
at least near the centre of the flow 22 = L/2 where we expect universal behaviour. 

In the preceding section we recognized another macroscopic length besides L entering 
the equation of motion: the marker trajectories s(t), for which there does not seem to be an 
analogue in critical phenomena. To be more precise, only the S2 component of s is relevant 
because of translational invariance in planes parallel to the plates. We refer to S2{t) as a 
"dynamic outer length" as opposed to L being a "geometric outer length". Because of its 
time dependence it cannot function as a correction scale in flip]), but it has to be dealt with 
during renormalization, cf. Sect. [TV]. 

We consider such an inverse RG to make more sense physically than a formal application 
of RG procedures for critical phenomena to turbulence. But this insight does not tell us 
how a procedure for eliminating degrees of freedom or a field-theoretic RG theory should be 
constructed. The main problem is our missing information about the statistical distribution 
for the r-velocity fluctuation. 



B. A new idea 

Instead of trying to construct a specific RG procedure based on some assumed statis- 
tical distribution, the new idea of this work is to make a general ansatz for possible RG 
transformations and to examine their scope under physically motivated constraints. 

Instead of specifying a procedure for eliminating degrees of freedom, we consider an "ar- 
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bitrary" transformation v(r, t) i— > v'(r', t') of the r-velocity fluctuation. But we demand that 
the renormalized r-velocity fluctuation v' obeys "essentially" the same equation of motion. 
This will reduce the number of transformations admitted but, in addition, we have to pose 
sufficiently many, physically motivated constraints which reduce the number of transforma- 
tions admitted even further such that nontrivial statements about scaling exponents become 
possible. Behind this idea is the hope that only a limited amount of RG transformations is 
compatible with the Navier-Stokes equation. 

The most important constraint to the equation of motion for v' is the invariance of the 
nonlinear term since the structure of the nonlinearity is prescribed by hydrodynamics and 
must not change. To give an example, because of Galilean invariance u-Vu must not be 
multiplied by a coupling parameter which changes under RG transformations. Although the 
r-velocity fluctuation v itself is Galilean invariant, such a variable coupling parameter would 
contradict the Navier-Stokes equation from which the equation of motion for v results. 

Since the markers' Eulerian velocity u(s, t) enters the equation of motion for v, the 
transformation of this quantity has to be specified, too. We demand that the equation of 
motion for u'(s',t') shall again be a Navier-Stokes equation which might differ from that 
for u only in the dissipative term. This is even more restricting than what we demand for 
the transformation of v. But we allow that v' does not have to be a difference of Eulerian 
velocities u' according to (|8|), because otherwise the transformation of v would be fully 
determined by that of u, and there would be no need and no gain to introduce v in the 
beginning. Since we believe velocity differences to be of fundamental importance, we look at 
the variables v and u as independent velocity fields which might "decouple" in the course 
of renormalization. Here, we do not see any analogy to critical phenomena. Note that the 
representation of the equation of motion (|TT| ) in terms of v, u, and U is unique if only 
Eulerian velocities at the current marker position s(t) enter. Although the reasoning for the 
invariance of the nonlinearity is based on the connection between the equation of motion for 
v and the Navier-Stokes equation, we keep this demand also if v and u decouple as outlined 
above. 

13 



The second big constraint comes from the assumed special role of energy dissipation. It 
seems natural to assign the central role which the free energy plays in RG theory of critical 



phenomena to the energy dissipation rate in turbulence p8| , |5T| . 

The energy dissipation rate (per mass) averaged over the full flow volume Q, 

£n = U 2 £5, Ir^Lir (^ (X) r) + ^ (Xj r) ) 2 (16) 

is a fundamental physical quantity in turbulence since it equals the energy injection rate 
£i n to maintain dynamical stationarity. We see it as the analogue to the free energy density 
in critical phenomena and expect it to be finite in the turbulent limit. This corresponds to 
the widespread belief that the limit — ^ is fundamentally different from the Eulerian case 
v = where Eq = 0. 

In RG theory of critical phenomena it is not the free energy density itself which is 
RG invariant but the total free energy. For turbulence, we suggest to also make such a 
distinction. Since e in = vUL(d X2 U(x 2 ))\ X2=0 / L (for plane shear flow) is determined by the 
mostly nonuniversal behaviour of the turbulent profile at the boundaries, does not seem to 
be well suited as a RG invariant quantity. Instead we choose the turbulent energy dissipation 
rate e t , which lies at the heart of Kolmogorov-based phenomenology, as an analogue to the 
free energy. To define it, we start from the mean local dissipation rate 

2 



e{x 2 ) = ^/(d x M i {x,t) + d Xi u j {x,t)') \ (17) 



near the centre of fluid motion (x 2 ~ L/2). The averaging procedure was defined in (^). 

We now rewrite e as a functional of r-velocity fluctuations and split off the dissipation 
related directly to the turbulent profile: e = E\j + e t with 

2 

(18) 



£v(x 2 ) = 7^(d X2 U(x 2 
£t(x 2 ) = - ( (d rj Vi{r,t) + d n Vj(r,t) 



1 

' r=0 

This splitting is unique since the mixed term involving both U and v vanishes. The profile's 
dissipation E\j is nonuniversal and corresponds to the background term present in the free 
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energy known from RG theory. But we do not know if the quantity J Q ^ £\j( x 2) is analytic 
in v as its correspondence to the background term in the free energy density, being analytic 
in the temperature, demands. We remark that e t equals the energy dissipation rate coming 
from the fluctuations u — U, i.e., 

e t (x 2 ) = ^({d Xj (u*(x,t) - Oi(x)) + d Xi ( Uj (x,t) - t/,(x))) 2 ^ . 

This turbulent energy dissipation rate e t shall be our "RG invariant quantity". This 
means two things, again analogous to RG theory for critical phenomena: first, any RG 
transformation should preserve the functional structure of e t when written as a functional of 
v(r, t) with parameters z/, x 2 , and L — except for an additive, v independent contribution. 
This corresponds to the functional invariance of the free energy density; again up to an addi- 
tive, order parameter independent contribution. The main difference to critical phenomena 
is the missing knowledge about the probability distribution functional for v(r, t). Instead, 
we work with the purely formal expression (|Tj|) for e t . Note that this might not be the same 
as using the definition (|TT| ) since v and u might be transformed differently. Second, and 
corresponding to the invariance of the free energy, the additive constant is put into e\j such 
that the renormalized value e' of the local dissipation rate equals the original value e: 

e' = e{j + e' t = eu + e t = e. (19) 

The aforementioned e* is now defined to be the turbulent limit of e t , i.e. e* = e t \ u ^o and 
is also assumed to be finite for v — > 0. For K41 scaling behaviour we identify the Kolmogorov 
length i] = (z/ 3 /^*) 1 / 4 with the dissipation length £, up to a dimensionless factor which is 
assumed to be a finite constant for v — > 0. If there are intermittency corrections, however, I 
should depend in a nontrivial way on z/, just as the correlation length scales nontrivially with 
the temperature parameter. On the other hand, the RG transformation behaviour should 
be £' = hi (analogous to a decrease of the critical correlation length by a factor b), whereas 
in general rj' ^ brj if v scales nontrivially. 

The other constraints are of minor importance and will be mentioned in the course of 
fully specifying our RG method in the following sections. 
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But before going into the details we can already infer some of the advantages and disad- 
vantages of such an RG approach. Since the forcing enters only structurally in the equation 
of motion for v, and we do not actually eliminate degrees of freedom, we cannot expect to 
be able to really calculate non-trivial scaling exponents. This is only achievable for critical 
phenomena since much more knowledge about the forcing, coming from thermodynamics, 
is made good use of. But we hope to learn if and how the lack of a boundary enforces the 
classical exponents and if there are universal relations between exponents. In our approach, 
we assume the validity of the "inverse RG" analogy to critical phenomena. The existence 
of RG transformations with the desired properties including consequences of symmetries 
is assumed, not shown by construction. Despite these disadvantages, we overcome several 
shortcomings of previous RG methods in turbulence. We do not have to know the proba- 
bility distribution for some velocity field or make wild assumptions about it. Our method 
is nonperturbative which is nice since fully developed turbulence has no obvious expansion 
parameter. 



IV. RG TRANSFORMATION 

We come to the actual presentation of our RG method for turbulence, following the spirit 
of the previous section. Because of the invariance of the nonlinearity the original equation 
of motion and its renormalized counterpart should have the following form (primes denote 
renormalized quantities): 

fyv + 5P[v-V r v] = £[v] = uS[A r v] + JF[v], (20) 

fyV + SVy-Vr'V] = C'[v'} = is'S[A r ,v'} + (21) 

Here, the terms linear in v, respectively v', are combined by introducing an inhomogeneous 
linear operator C, respectively Renormalization of forcing and dissipation expresses itself 
as £ 7^ C Besides v 1 > u, one expects to find structurally new terms, which we initially 
consider as part of T' ^ T . They have to be interpreted as 'dissipative' or 'driving' on an 
individual basis. 
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In order for equation (pTj) to make sense, the transformation of r, t, v, s(t), and u 
has to be defined. Finding transformations such that fl2~T|) has the prescribed structure is 
the main goal of Sects. [TV] and [V]. We begin with discussing the parameters which these 
transformations may depend on. 

An important assumption is the regularity of RG transformations for v — > 0. For critical 
phenomena, the analogous regularity in the temperature approaching the critical temper- 
ature is a main achievement of RG theory. In this work we concentrate on the universal 
scaling behaviour in the inertial range and thus may already carry out the limit v — > in 
the transformation. Effectively, it suffices to take v independent transformations. As long 
as we do not look at multiple-time correlations of v, we assume that the transformations 
do not have to be t dependent. But r should be admitted, as well as both outer lengths L 
and S2(t) since all these are nonuniversal. Dismissing s± and S3 takes into account geometric 
symmetries, and the possible dependence on S2(t) allows an implicit time dependence. We 
even start with the admission of all components of s in order to study the consequences of 
the restriction to S2 later. Note, however, that we will not enforce this restriction for the 
transformation of u; it will be allowed to depend on all components of s. With the prelude 
in mind, we come to the specification of our ansatz for RG transformations. 

A. Space and time 

Spatial distances are just rescaled: 

r' = (1 + a)r + 0{a 2 ), thus d r > = (1 - a)d n + 0{a 2 ). (22) 

The RG rescaling factor b = 1 + a deviates only infinitesimally from 1 by the RG flow 
parameter a. Experience with critical phenomena teaches that such infinitesimal RG trans- 
formations usually exist and are even easier to deal with. Because the RG transformation 
shall be inverse, we have r' > r. For time differences we will also begin with a simple rescaling 
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ansatz, 



t' = t + a(t-t )z + O(a 2 ), 



(23) 



with a single constant scaling exponent z. In a supplement to Sect. |V A| we will examine a 
more general transformation of r and spatially varying z. 



B. Marker trajectories 

Because of the presence of the marker trajectories s(t|x , to) m the equation of motion we 
have to specify how they shall transform in spite of not knowing a corresponding quantity in 
critical phenomena. Since the trajectories are nonuniversal their transformation could and 
should be different from the transformation (22) of distances r. 



Because the trajectory's position s(t) enters the argument of the Green function in (12), 



its transformation is linked to that of the geometry. We demand that the geometry, and thus 
the Green function, is RG invariant. This corresponds to keeping the geometric structure of 
the microscopic lattice fixed in critical phenomena. On the other hand, the RG transforma- 
tion should change the nature of the forcing as it changes the microscopic interaction for 
critical phenomena. Thus, keeping the geometry invariant actually assumes that all renor- 
malization of the forcing can be covered by a change of the linear forcing T . Because of this 
invariance of geometry and V = L we have to define V in ([21]) to be 

£>Mr, t) = ft (r, t) - f d 3 p' d r ,G n (r' + s'(t'), p' + s'(t')) ^(p, t). (24) 

Jn(s'(t')) 

The most simple way to satisfy this geometric invariance is to keep the marker's position 
the same, i.e., s'(t') = s(t). For a fixed fluid element serving as the marker this is impossible 
since the dynamics and the trajectory of a fixed marker change inevitably under RG trans- 
formations. But we may allow the marking fluid element to change. We choose, for every 
time t, a marker with initial position x such that 

S(t) EE S (t|xo, t ) = S(t'|x , t ) EE S \t'). (25) 
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This definition is unique because of the one-to-one correspondence of initial positions xo and 
current positions s(t|x ,io) f° r ai l t. It also guarantees s'(t') G f2. The set of all transformed 
trajectories s' naturally defines an Eulerian velocity field 



u'(s',0 =s> = d t ,s'(t'\x' ,t )L 



(26) 



which should satisfy a Navier-Stokes equation. We will cover the transformation of u in 



Sect. VI B 



Going from xo to Xq enforces to also relate the transformed fluctuation v' to a marker 
with initial position Xq instead of xq. Therefore the time derivative dfV in ([H]) has to be 



taken with fixed Xq. From fl23|) , (|25 ) and ( 26 ) we then find 



ir ,x 



:i-az)d t \ sr + u'.V s \ rt + 0(a 2 ). 



(27) 



C. Statistical averages 

Mean values (/) are defined analogously to (I7|) with the corresponding transformed 
quantities. For the transformations defined so far we simply have (/) = (/). To see this, we 
only have to look at the spatial averaging: 



dh 



/3 



2\s' l (r')-x l \<a' « 

(i 3 x 

2\si{r)-x i \<a ^ 



dh 



-0 



2\s t (T)-Xi\<a' & 

1 _ 3a + a ^J2 [<%-^o ; ;-f) + ^i-aroy+f)] + 0{o* 



Here we made use of (|25l) and a' = (1 + a)a + 0(a 2 ). The latter holds because a is a 
microscopic length which should be rescaled within an inverse RG theory with a factor 
b (with no nontrivial exponent, cf. Sect. |III A| ) - - independent of possible intermittency 
corrections because we identify this formal parameter a with a physical quantity, namely 
the linear size of fluid elements, which we also designate as the hydrodynamic length scale. 
The formal limit a —>■ corresponds to a -C £, i.e., one is only interested in the hydrodynamic 
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length being smaller than the dissipative scale. This role of a is analogous to that of the 
system size L in critical phenomena, where the case L exceeding the correlation length is 
of main interest. Just as finite size effects do not change the critical exponents, we expect 
intermittency exponents to be independent of effects on scales less than a. 

We now assume that the surface averages represented by the terms with a delta function 
equal the full average (), at least after taking the time average and limits according to (ff|). 
Then all terms of order O(a) cancel and we arrive at (/) = (/) + 0(a 2 ). 



D. r-velocity fluctuations 

The crucial part of our ansatz is the transformation for the r-velocity fluctuations map- 
ping v(r, t\s) to v'(r', t'\s'). This shall be realized through a linear transformation operator T 
which expresses the abstract elimination of degrees of freedom and the subsequent rescaling: 



v' = v + aT[v] + 0(a 2 ). 



(28) 



The linearity of T avoids terms more-than-quadratic in v in the equation of motion. T will 
in general not be a simple rescaling, but may transform v and its derivatives differently. 
This takes care of the possibility that different orders of turbulent velocity derivatives might 
be independent scaling fields. Such additional scaling fields would be analogous to the addi- 
tional order parameters in dynamical critical phenomena introduced by conserved quantities 
relevant on the large scales. 

Of course, we have to specify an explicit expression for T. It shall cover all linear differ- 
ential operators of finite order in r, s and L. This choice is not motivated physically. 

JV M 

L,( n > m ) Pi . . . PI Pim _l a 

i,j;ki,... ,k n s k t 



UL i,jM,-,k n u r kl 



f) f) m _l d n ' m ) Pi . ..Pi Pfm 



sm = E E 

n=0 m=0 

The derivatives d Si and di are defined with r kept fixed according to 



Vj . 



(29) 



V S = V S 



r,t,L 



dj, = d T 



I r,t,s 



:v xo s| , L r i -v X0 \ 



d Ar,t,x - d L S \t,x - V s\ r ,t,L- 



(30) 
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This is well-defined since the uniqueness of the marker trajectories allows to unambiguously 
rewrite v as a function of s instead of x . It should be stressed that the time derivative in 
the equation of motion for v and w is not taken with constant s, but xo, i.e. dt = oU r . 
The characteristic property of v (and w) being velocity differences, cf. allows to express 
s derivatives through r derivatives: 

<9 Si v = S[d r .v], d Si d rj v = d r .d rj v, etc. (31) 

Because of this "difference property" mixed derivatives in r and s are not independent 
and thus omitted from (p9|) . The coefficients a^™'.^ kn and ftjj'™^ kn are (continuously 
differentiable) functions of r, L, and s (later s 2 only)- Without loss of generality we can 
set f3?J m) = and assume these coefficients to be symmetric in the indices k\, . . . ,k n . The 
finiteness of M and N will not be used explicitly, but in order to find or exclude solutions for 
T with infinitely many contributions in (p9[) , we would have to study questions of convergence 
and exchange of limits in a suitably chosen function space. 

It seems reasonable that v' should have the same basic properties as v. Obviously v' has 
to be divergence free and has to vanish for r' = 0. We also demand a zero mean of v' as a 
consequence of homogeneity. Furthermore, we try to preserve the difference property. But if 



u' and v' decouple as argued in Sect. |III B| this constraint could or should be relaxed. These 



four properties lead to corresponding constraints for T: 

V r -T[v] = for V r -v = (divergence freeness), (32) 

T[v]| r=0 = for v(0,*) = (origin), (33) 

(T[v]) = for (v) = (homogeneity), (34) 

(SV r - V s )T[v] = SV r v (35) 
for V s v = ^[Vrv] (difference characteristic). 

They are not automatically fulfilled by ansatz (^) and will be enforced later (Sect. |V A| ). 
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E. Eulerian velocities 



Since the Eulerian velocity u' defined in (^6]) should follow a Navier-Stokes equation and 
because of the formal similarity to the equation of motion for v, we employ a transformation 
T for u analogous to T: 

u' = u + af[u] + 0{a 2 ), 

N M (36) 
n=0 m=0 

The coefficients dtf*'™^ kn are functions of s and L and each is symmetric in the indices 
ki, . . . ,k n without loss of generality. Divergence freeness of u' leads to the constraint 

V s -T[u]=0 for V s -u = 0. (37) 



V. INVARIANCE OF THE NONLINEARITY 

After having specified the transformation, we can apply it to the equation of motion for 
v. But at first, only the linearity of T is used, and we introduce the abbreviations V = ST> 
and V = SV. After inserting (0), fl27|) and fl2|) into the lhs. of (0) and d t \ from @ we 
find 

dfV + P'[v'- V r /v'] = (l-az + aT)C[\] + aT[v] + (u - (1 - az)u)- V s v 

+ a(z- 1 -T)P[vV r v] +a5V[v-V r v] + aV[T[v]-V r v + v-V r T[v]] + £>((t 2 ) 

^r'fv'i + o^ 2 ). 

Here, is defined through V' = V + a5V + 0(a 2 ) and T in an abbreviation for d{T — Tdt- 
The nonlinear terms have to cancel out, since CJ is supposed to be inhomogeneous linear in 
v: 

P[T[v]-V r v + vV r T[v]] -T?[vV r v] =-(( 2 -l)? + 5P)[vV r v]. (38) 
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This inhomogeneous linear constraint for the operator T is at the heart of our RG approach 
for turbulence. We call such a constraint the "nonlinearity constraint" . From the remaining 
terms and because of C'[v'] = C'[v] +o~£T[v] + 0(a 2 ) we infer an equation which determines 

£'[v] = (1 - az)C[w] + a(TC - £T)[v] + af[v] + (u' - (1 - az)u)- V s v + 0(a 2 ). (39) 



It will be analyzed in Sect. [VJ. The inhomogeneous linearity of £ will then be confirmed a 
posteriori. 

We were not able to find the general solution of the nonlinearity constraint (|38|) within 
the ansatz fl2~9"|). The main difficulties arise from the nonlocality of the projector V. We 
therefore first look for solutions which fulfill a simpler constraint, namely we ignore the 
pressure term. The resulting modification of the nonlinearity constraint allows us to find its 
general solution. Formally this means to substitute the identity operator for V and to set 
5V to zero in (RSI) . This leads to the considerably simpler constraint 



T[v]-V r v + v-V r T[v] -T[v-V r v] = -(z - l)v-V r v. (40) 
In a second step, we demand that the solution T of (HO) meets the full nonlinearity constraint 



(pg). Thus, the effects of pressure are only taken care of subsequently. For this second step 
it is sufficient that 

(TV - VT)[q\ = 5V[q\ for q = v-V r v. (41) 

In general, this constraint will reduce the set of solutions. The geometry of the flow volume 
enters only here. One should keep in mind that this two-stage method might not find all 
solutions of (p8|). For example, there might exist additional solutions where the pressure 
term and the local nonlinear term "mix" or the equality in fl4*0| ) holds only except for a 
divergence free, additive term. 
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A. Neglecting pressure 



Determining the solutions of ( HDp is quite tedious and boring. The lhs. is calculated for 
all terms of the sum in ansatz fl29l) . It is assumed that all components of v and all their 
derivatives are independent from each other — except for derivatives which are coupled 
because of the divergence freeness and must be reexpressed through independent derivatives. 
Then, fl4~0|) must be individually fulfilled by all terms with a common derivative structure. 
More calculational details of this straightforward procedure can be found in the Appendix 
[X[ Here, we proceed to the result: The general solution of (f40|) may be written in the form 

%[v] = oiij vj - [(z - l)rj + a jik r k - ^ij\{d rj Vi) + Pjd Sj Vi + (3 L d L Vi (42) 

with 17 free parameters a^j, jj, j3j, /3l, and z, which are arbitrarily depending on s and L 
(except for z = const.). The a^j are dimensionless, whereas 7^, (3j and (3l carry the dimension 
of a length. It should be noted that divergence freeness is not relevant for this result: The 
same general solution holds if one assumes all derivatives to be independent of each other. 

Now we have to check the constraints (|32"l) -(j35|). The divergence constraint ([32]) and 
the homogeneity constraint ([34]) are automatically fulfilled by all transformations (f42D - The 
origin constraint (|33|) demands jj = for all j = 1, 2, 3. The difference constraint (|35D leads 
to 

(Sd. n - d Sl )%{\] = -[(z - l)5ji + a jt i + (<9 Si /3 j )]5 , a r .w i + Ijdr^v^^ 

- {{d Sl Ui,j)vj - [d Sl ((z - l)5 jk + a jtk ))r k d r .Vi + (d sl jj)d r .Vi + {d Sl (3 L )d L Vi} 

= Sd ri Vi. 

The terms in curly brackets all have a different structure and must vanish individually. 
Consequently, all coefficient functions except (3j must be independent of s, i.e., they are 
free constant numbers, which is compatible with the constantness of z. From the remaining 
terms we find again jj = and a linear dependence of j3j on s: 

Pj = —z Sj — ctjj si + Po-j with p .j = const. (43) 
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We see that the difference constraint restricts the admitted transformations considerably. 
This goes as far as enforcing a dependence on all components of s which contradicts our 
desire of a sole S2 dependence. Later on, we will therefore also consider weaker forms of the 
difference constraint. In the extreme case of ignoring this constraint we have to generalize 
the ansatz ( |2Tf ) to include mixed derivatives with respect to s and r. But this does not lead 
to a more general solution (fl2"|); the difference constraint does not change the solutions for 
T in the pressureless case. 

Supplement: anisotropic rescaling Imagining a definite RG transformation with elimi- 
nation of degrees of freedom in a thin boundary layer, an anisotropic rescaling of spatial 
distances seems much more appropriate than the isotropic rescaling adopted from critical 
phenomena. We therefore repeat the preceding calculations with 

r ■ = n + o-RijTj + 0(a 2 ), R iyj = const., 

instead of (p2|). Since a is only determined up to an arbitrary constant factor, we may choose 
to fix one of the parameters Rij, e.g. i?2,2 — 1- This general ansatz includes the physically 
motivated case R it j = ^,2^,2- We indeed find a solution generalizing (|X5]): 

af 1 = — (ocij ^ — Rj,i)rj — zri + with 7; independent of r. 

But the full solution for T corresponding to (^) complies with the divergence constraint 
(j32|) only for Ri j = 8 it j. We conclude that the isotropic nature of incompressibility enforces 
an isotropic rescaling of r, at least for the set of transformations found so far. 

Supplement: spatially varying z exponent A crucial but perhaps too restrictive assump- 
tion is the necessity of only a single dynamical exponent z, which might be a speciality of 
critical phenomena and not valid for turbulence. We therefore might try to use a position 
dependent exponent z(s(t)/L) exploiting the fact that there seems to be no quantity anal- 
ogous to s in critical phenomena. Without elaborating on all details (cf. Ref. here we 
just mention that such a generalized z is compatible with the transformations T considered 
so far, but implies homogeneity and difference constraints which cannot be trivially satisfied. 
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The main problem, though, comes from the nonlinearity constraint upon the transformation 
T for u which will be dealt with in Sect. |V1B| . It enforces z = const., and so we shall not 
pursue this idea further. 



B. Invariance of the pressure term 

Now we perform the second step where we have to make sure the commutator constraint 
(HI) holds. One verifies very easily that T commutes with S (since jj = 0), and it suffices 
to check to exchangeability of T with T> under the action of S: 

S(TV - VT - 8V) [q] = for q = v- V r v. (44) 

Here, we define 8T> through T>' = T> + cr 5T> + 0(a 2 ). After integration by parts we find from 



([jjD and fl42|) for any Green's function Gq 

(TV-VT- 5V)i[q\ = - [ d 3 p { (a itj + a jA ) (d rj G Q (r + s, p + s)) 

Jn(s) L 

- d n (z {r j d Tj + pjdp. + 1) + a jtk (r k d r . + p k d Pj + S jk ) 

- /3jd s . - (3 L di)jGn{j + s, p + s)j | (d pi qi(p,t)). (45) 

This expression can now be evaluated for different Green's functions corresponding to differ- 
ent geometries. Although we would like to focus on plane shear flow because of its physical 
significance, we will also study an "infinite" system which is similar to what is considered 
in most other theories of turbulence. This strategy enables us to get insight about the con- 
sequences of the presence of boundaries. 

Let us first consider the infinite system's Green's Function, 

Gn(x, » G f (x, = -^j^jy (46) 

We prefer to call this substitution an "approximation" to the plane shear case since a 
truly infinite turbulent system does not make sense within our framework with its physical 
forcing. Effects due to the boundary conditions and the finiteness in one space direction are 
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neglected. Therefore, we keep the transformation ( p5|) of s for the infinite system although 
the reasoning for it is not transferable since T> is independent of s. With the help of the 
relation 

(r j d rj +p j d Pj + l)G i (r,p) = (47) 
the commutator (f|5|) can be evaluated as 

(r k ~ Pk)(rj - Pj) 



(TV -VT- 5V)M = - [ d 3 p (3 VW + a kJ ) ^—^ 



2 



- «iJ (l " 3 } (^A(r, p)) (0„ ft (p, f)). 

The commutator constraint (H) holds only if 

a^fc = — otkj for j ^ k and independent of j. (48) 

This expresses the rotational symmetry of the problem. 

Translational invariance and the absence of a length scale set by the boundary conditions 
are responsible for the independence of G{(r + s, p + s) from s and L. We demand this 
independence also for the transformation T: 

azij,z independent of s, L and (3j = /3l = 0. 

This contradicts the difference constraint ( f4"3|) (except for the uninteresting case a^j = 
z = 0). In order to avoid s dependent transformations, one has to give up the difference 
constraint. This seems reasonable with respect to the decoupling of v and u discussed in 
Sect. [Ill BL Then, the general solution fl4"2]) may be written in the form 



7~[v] = ctij Vj - axTjd^Vi - - ^ a j>k (r k d rj - rjd rk )vi, (49) 
ai = z — 1 + ao, o>q = ctjj, (50) 
with 5 independent, constant and dimensionless parameters 0-1,2, «i,3, 0:2,3, a i, an d z. 
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Now we turn to the physically relevant case of plane shear flow. For its Green's function 
Gn, Eq. ©, 

( r A + Pj d Pj + s 2d S2 + Ld L + l)G n (r + s, p + s) = (51) 

holds instead of (^7|). The complete evaluation of the commutator constraint ( ^5| ) leads to a 
rather complicated expression for which it is difficult to find all solutions. But we find (f45"D 
to be satisfied by all transformations T with 

ol\,2 = a 2 ,i = a 2 ,3 = c*3,2 = 0, ai >3 = -013,1, 

(52) 

a 1 as in (|50|), /3 2 = — (1 + ai)s 2 , /?l = -(1 + «i)L. 

For this, one only needs to know that Gq has to adopt the form Gq(t + s, p+s) = Gq((t+ 
s)/L, (p+s)/L) jh because of dimensional reasons and that Gq depends on n, pi, r 3 , and p 3 
only in the combination (r x — pi) 2 + (r 3 — p 3 ) 2 because of translational invariance in planes 
parallel to the plates. 

Due to this symmetry we additionally demand that the transformation should not depend 
on s\ and s 3 as discussed before. This allows for explicitly space dependent transformations 
with 

aij(s 2 /L) independent of Si,s 3 and /3i = /3 3 = 0. 

But this space dependence again contradicts the difference constraint. Once more, we decide 
to give up this constraint in favour of obtaining a different set of admissed transformations 
which can be written in the form 

%\y] =a Vi + ai l3 (5j ) itf 3 - 5i. 3 t>i) 

(53) 

- axrjd^Vi - cti i3 (r 3 <9 ri - r x d r3 )vi - (1 + «i)(s 2 <9 S2 + Ld L )vi 

with 3 independent and dimensionless parameters, namely qji i3 and oti, depending on s 2 /L, 
and constant z. Unfortunately, we can not come up with a satisfactory physical interpretation 
of this transformation yet. 
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VI. TRANSFORMATION OF DISSIPATION AND FORCING 



Having found sets of admitted transformations, (^9|) and (|53|) respectively, we transform 
the linear part £' of the equation of motion for v. 

A. Viscosity 

We begin with the calculation of (1 — az)£ + a{T £ — £T) in ( 153 ) for the viscous part 
vSA r of £ in the case of plane shear flow: 

(1 - (Tz)vSA r v + a(TvSA r - uSA r T)w = (1 + a(2a 1 - z))vSA r v + 0(a 2 ). 

Taking into account the rescaling of A r we find u'A r >V to be a renormalized viscous con- 
tribution to C'[v'] with the rescaled viscosity 

z/ = (l + („cr)is + 0(a 2 ) with (u = 2(1 + a x ) - z. (54) 

Here, the symmetry properties of the parameters aij, cf. (^8|) and (0), have to be used 
explicitly. This result also holds for z(s2/L) and also for the infinite system approximation. 
In the former case v has to be generalized to a function of S2/L. 

Since the viscous part of £ is reproduced with exactly the same structure, we can safely 
regard all other terms in (^) as contributions to the renormalized forcing T' . The property 
v' > v of an inverse RG transformation implies the constraint ( u > 0. 

B. Forcing 

The admitted transformations T for the Eulerian velocity u which are covered by the 
ansatz (|36|) can be deduced easily from the results of Sect. because of the formal similarity 
of the Navier-Stokes equation and the equation of motion for v: 

%[u) = a itj Uj - [z sj + a j)k s k ](d Sj Ui) + (3 L d L Ui. (55) 
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The parameters atj± and $l are geometry dependent in the same way as a^j and (3l, but in 
any case have to be independent of s. In general, the viscosity in the Navier-Stokes equation 
rescales with a different scaling factor and therefore will be denoted with the new symbol v: 

v' = (1 + ( v a)v + 0(a 2 ) with ( v = 2(1 + 5i) - z. (56) 

Next, we calculate the last two terms in (^). First we look at T for plane shear flow. In 
the most simple case of constant parameters ao, oti, and we find from (|53|) 



t = d t T- Td t = -(1 + ai) [u 2 d S2 - s 2 (d S2 u)-V s - L(d L u)-V s ] . (57) 
Here, we made use of the commutators <9 t V s — V s <9i = — ( V s u) • V s and d t di — did t = 



(9lu)-V s which follow from definition (]30|) . We see that aT has the same structure as 



(u ; — (1 — az)u)- V s in (|39f) , but these two contributions to T' have a different structure 
than the remaining contribution <r(TjF — JFT), since the former do not contain the turbulent 
profile U. Thus, they describe structurally new, additional forces (per mass). Since such an 
additional forcing is compatible with the idea of our RG approach, one should study how it 
behaves under the RG transformations' next iteration. We now show for one specific example 
that structurally new terms appear with each RG iteration. A term in ([57]) like S2(9 S2 u)-V s 
in T results in a contribution 



-a(l + a 1 )(s 2 ) 2 [(d S2 d S2 u)-V s -(d S2 u)-V s d. 



S 2 J 



to T' from <t(TJF — JFT). One easily sees that with each iteration a new contribution with 
an additional s 2 derivative is created. The same holds true for nonconstant a parameters 
which make appear even more distinct forcing terms. One might try to somehow get rid of 
this additional forcing. We thoroughly analyzed several possibilities which all resulted in the 
necessity to give up one or the other fundamental requirement pPf . We conclude that the 



appearance of additional forcing terms is inevitable, although we could not further analyze 
the transformation behaviour. 



30 



We are now able to calculate the turbulent profile, that is the mean value of the trans- 
formation ( |55"D for u. For plane shear flow we have U' 2 = and 

U' l = U i + cr{a^U j -(ao + z)x 2 d X2 U i + p L d L U i } + 0(a 2 ), (i = 1,3). (58) 

Unfortunately, this equation does not allow a unique determination of the profile. But all 
profiles which take on the functional form U* = Uif{x 2 /L) and are bounded for all x 2 and 
v — > are admissible solutions of (|58|) (with «o = 0). In the remaining part of this work, we 
will consider only this case «o = 0. 

Finally we are in the position to discuss the total transformation of the forcing T. The 
calculation of the remaining contribution (1 — az)T + o{TT — FT) in (|39|) shows that T 
would be structurally reproduced (except for the additional forcing discussed above), if 

U'i = Ui + cr{a itJ Uj - (« + z)x 2 d X2 U l + /3 L d L Ui} + 0{a 2 ) (i = 1, 3) 

would hold. Since this is not true (cf. (|58"D) yet another, U dependent type of additional 
forces appears. 

We conclude that a closed, iterable RG equation for the forcing T does not seem to exist. 
Clearly, more work on the transformation properties of the forcing is necessary. For example, 
one would have to examine if the additional forces are irrelevant in the sense of RG theory. 
Also, a restriction of U such that the iteration becomes closed could be possible. 

VII. ENERGY DISSIPATION AND STRUCTURE FUNCTIONS 

We are done with the transformation of the equations of motion for v and u. This 
enables us to look how physical quantities, more specifically energy dissipation and structure 
functions, behave under the set of admitted RG transformations for v, namely (|49| ) with the 
parameters a\ t 2, ai,3, 02,3, Qti, and z for the infinite system and ([53]) with the parameters 
a li3 , cti, and z for plane shear flow, respectively. 
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A. Invariance of energy dissipation 



The transformed turbulent dissipation rate e t is obtained from the representation 
and the RG transformation (K), depending on the choice for the parameters atij, j3j, and 
Pl. 

For the infinite system approximation our requirement of structural invariance of e t is 
fulfilled automatically, except for a simple rescaling: 

e[{y') = (l + ( E cr)e t (u) + 0(a 2 ) with ( e = £„ + 2(a - 1 - a x ). (59) 

Structural invariance is most naturally satisfied by setting ( £ = 0, and this also keeps the 
value of e t the same: 

4 = 4K) =£t(V) =e t [y) =£ t . 

Because of our second invariance constraint for the energy dissipation, e' = e, the value of 
the profile's dissipation is invariant as well, e[j = E\j. But we do not have any knowledge 
about the functional structure of s'^v') as long as we do not know how to handle the 
transformation of the forcing. 

For plane shear flow the structure of e t does not reproduce itself: 

e' t (u', x 2 , L, U') = [1 + Ce<r - a(l + ai )(x 2 d X2 + Ld L )\ s t (u, x 2 , L, U) + 0(a 2 ). (60) 



The contributions in addition to ( |5*9]) are interpreted as nonuniversal since they are related 
to the boundaries of the system, and are therefore put into e' v : 

e' v = e' v (v', x 2 , L, U') + a(l + a 1 )(x 2 d X2 + Ld L ) e t (v, x 2 , L) + 0(a 2 ). 

Again we set ( £ = to guarantee structural invariance of e t , i.e., e t {v') = St(v). 

In both cases the condition ( £ = serves as an additional constraint for the RG trans- 
formations T which allows to eliminate one parameter which we choose to be a$. Using 
( u from (|54]) and Oq from (|50|) we find Oq = |(1 + oti) and z — |(1 + «i) with the only 



remaining free parameter a\. This value for z gives us C, v — g(l + ai), which is different 
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from ( u = |(1 + cx,x). This means that u, the kinematic viscosity damping the r- velocity 
fluctuations, and P, responsible for viscous losses of the Euler field, scale differently. 

B. Scaling behaviour of structure functions 

Finally we come to the transformation of the structure functions D^ m \ defined in (|HD, 
and their inertial range scaling. Contrary to e t we demand the same functional form for 
plane shear flow and the infinite system, i.e., D^ m ''(r',iy',X2,L) = D^ m '(r',p',X2,L), thus 
not allowing for an additive renormalization. 

Again we start with the infinite system and find for the set of admitted transformations 

£M(r', v>) = {1 + a[ma - (a^ + J2^, k r k )d rj ]}D^\v, v) + 0(a 2 ). (61) 

We seek for scaling solutions to this RG equation in the limit v — > which obey the geometric 
symmetries of the infinite system, i.e., depend on r = |r| only. Then, ([Jl]) simplifies to 

D (m) (r', v') = {1 + a[ma - aird r )}D^ m \r, v) + O[o 2 ). 

Because L is assumed to be the only scale for intermittency corrections, only solutions are 
sought which are bounded for v — > with r ^ fixed. Using r' = (1 + cr)r on the lhs., the 
unique solution then is 

D^rJocrC- with ( m = = ™ (62 ) 

1 + ax 6 

It is not surprising that K41 behaviour is enforced since the admitted RG transformations 
for v turned out to be independent of L. But this result at least confirms consistency of all 
assumptions behind our method. 

For plane shear flow the RG equation 

D ( - m \r',i/,x 2 ,L) = {l+a[ma - (airj + ^ a i,k r k)d r] 

- (1 + ai)(x 2 d X2 + Ld L )]}D^ m \r, u, x 2 , L) + 0(a 2 ) (63) 
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holds. Here, respecting again the geometric symmetries, we are looking for solutions depend- 
ing on r, x 2 , and L alone, thus simplifying flB3| ) to 

D { - m \r',v',x 2 ,L) = {1 + a[ma - ct x rd T 

- (1 + ai )(x 2 d X2 + Ld L )]}D (m \r, v, x 2 , L) + 0(a 2 ). 

All solutions bounded for v — > can be represented in this limit as 

D^ m \r,x 2 ,L) oc r m/3 f c (r/L,x 2 /L), with arbitrary function / c . (64) 

This allows for scaling solutions of the type _D( m ) oc r^ m with arbitrary inertial range ex- 
ponents ( m . Within our method we could not yet identify more constraints which would 
restrict the exponents or even nail them down to a specific value. 

VIII. SUMMARY AND DISCUSSION 

A. Summary of the calculations 

The results of Sects. [TV] to |V11| can be summarized as follows: with a quite general ansatz 
for the transformation T for the r-velocity fluctuation v we evaluate the nonlinearity con- 
straint originating in the invariance of the nonlinear term of the equation of motion. First 
this is done ignoring the pressure term which restricts T to a class of solutions with 17 
parameters which are allowed to depend on s(t) and L (except for the dynamical exponent z 
which must be constant). All further conclusions were drawn for this set of transformations. 
The difference constraint which expresses that v is a difference of two Eulerian velocities 
is abandoned since it turns out to be incompatible with a dependence of T on s(t) accord- 
ing to geometric symmetries. This is somewhat surprising and means that the renormalized 
r-velocity fluctuations decouple from the Euler field to become quantities of their own signif- 
icance. Taking into account a commutator constraint for T which guarantees the structural 
invariance of the pressure term and some other constraints restricts the set of solutions fur- 
ther. The number of free parameters is reduced to five constants ai >2 , 0^3, 02,3, «i, and z 
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for the infinite system and to one constant z and two functions and a,\ of S2/L for plane 
shear flow. 

The invariance of the geometry under renormalization suggests a simple transformation of 
the marker trajectories which may be loosely referred to as s' = s. We find a transformation 
T for the Eulerian velocity u which preserves the structure of the nonlinearity of the Navier- 
Stokes equation. It introduces four new constant parameters for the infinite system and two 
new constant parameters for plane shear flow. The RG equation for the turbulent shear 
profile advises a special choice for T. The viscous term of both equations of motion for v 
and u reproduces structurally, but with differently rescaled viscosities v respectively v. But 
for the forcing T several types of structurally new terms appear whose physical meaning 
is not yet understood. We accept such additional forcing, but do not examine it further at 
present. Neither the specification of a closed RG equation for T nor even its existence can 
be established as matters stand, and we only consider solutions T with constant parameters. 
The second important restriction coming from the invariance of the energy dissipation rate 
is fulfilled by setting the rescaling exponent of the structurally reproducing part of the 
turbulent dissipation rate to zero. This eliminates one parameter of T. 

The RG equations for the structure functions are examined for inertial range solutions 
with prescribed symmetries corresponding to the geometry of the flow. For the infinite system 
only K41 scaling exponents turns out to be allowed, whereas for plane shear flow there are 
scaling solutions with yet arbitrary exponents. This is an indication that intermittency 
corrections really require the existence of an outer length scale. But neither can the scaling 
exponents be calculated with this RG method, at least not in its present stage, nor can 
a physical mechanism for inertial range intermittency be identified. We also investigated 
periodic boundary conditions and half-space geometry |JIJ and found that again non-K41 
exponents are allowed, implying that both either a zero v-boundary or a finite distance 
between the plates are sufficient for the possibility of intermittency corrections. 
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B. Discussion of the method 



A main advantage of our RG method is its use of the Navier-Stokes equation and the 
question arises which of its aspects are essential for our results. 

The different physical meaning of Eulerian velocities and r-velocity fluctuations can be 
seen in the difference of the transformations T and T. The appearance of s in Green's 
function motivates the choice of the transformation of s and helps determining fa. But the 
occurrence of the operator S, which introduces the r-increment, cf. (0), in the equation 
of motion for v effectively does not enter the evaluation of the commutator constraint. 
Altogether one sees that it is not essential to introduce Monin's transformation for the 
results obtained so far. 

The quadratic and first-order gradient nature of the nonlinear term influences directly the 
admissibility of transformations through the nonlinearity constraint. The tensor structure 
and the dimensionality have no qualitative input on the RG transformation. Formally, all 
results also apply for general d dimensional turbulence, but since some assumptions such 
as the invariance of energy dissipation must not be valid in d = 1 and d = 2 dimensions 
the physical significance of the results is questionable. We stick to d = 3 in this work. The 
simplified approach of Sect. V A| to solve the nonlinearity constraint and the local ansatz for 
T do not take into account the strong nonlocality of pressure. It is one of the consequences 
of Galilean invariance that the fixed prefactor of the nonlinearity is responsible for the 
elimination of one of the transformation's parameters, namely ot\. Otherwise the scaling 
exponents would be arbitrary also for the infinite system. Concerning Green's function for 
plane shear flow, only its symmetry properties were utilized. We consider this as a signature 
of universality: the results are independent of the detailed nature of the boundary conditions. 

Incompressibility does not change the set of solutions found for T (Sect. [V A|) , but it 
is exploited in the definition of statistical averaging and enforces an isotropic rescaling of 
r (Sect. [V A| ). The definition of statistical averaging is purely formal, but the microscopic 
nature of the space average enters. The transformations found fulfill (v) = 0, but do not 
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take any special advantage of this. 

Now we mention some of the open questions which could be investigated further within 
this method. Most importantly, the transformation of the forcing should be worked out fully 
by finding a suitable treatment for the structurally new contributions. 

One could also think of looking for strongly nonlocal transformations not included in 
the ansatz (f?9|). We believe this to be quite difficult within the present framework since the 
simplified two-stage process for solving the nonlinearity constraint ( p8|) is not reasonable 
then. But if this is accomplished, it would again allow to consider an anisotropic spatial 
rescaling, new transformations for s, and an s-dependent dynamical exponent z, to give just 
a few examples. We cannot rule out that new transformations exhibit non-K41 scaling even 
for the infinite system, although this does not seem very likely. 

Of course, the transformation properties of other physical quantities could be studied. 
Top candidates are the dissipation length I and the energy dissipation correlation function. 
In principle, the examination of the scaling behaviour for v > could provide insight to the 
crossover from the inertial to the dissipative range. Finding v dependent transformations 
should also confirm our assumption of v independence in the turbulent limit. 

The application of our method to other dynamical systems whose scaling behaviour is 
better known, such as Burgers' equation, Kraichnan's turbulence model for passive scalars, 
and the GOY cascade model could be possible, but is not self-evident, since the physical 
idea and the key invariance assumptions do not seem to be transferable in a meaningful way. 

To summarize: although we could not calculate concrete numbers for intermittency ex- 
ponents, we found apparently non-trivial RG transformations which relate intermittency 
corrections to the existence of an outer length scale. Several open questions remain which 
should be investigated in the future. 

Acknowledgments: We thank Detlef Lohse for stimulating discussions and Stefan 
Thomae for bringing Monin's transformation to our attention. 
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APPENDIX A: SOLVING THE NONLINEARITY CONSTRAINT 



Here we present more details how to determine the RG transformations found in 
Sect. FSl 



The ansatz (^3) is plugged into the lhs. of the weakened nonlinearity constraint (TO) and 
derivations are performed such that all derivatives operate directly on v components only. 
This yields the following contributions to the lhs. of (fHf ): 

coefficient functions oifj : 

a ?f v A d nVi) + {d^ff^Vj (Al) 

coefficient functions : 

coefficient functions cufj : 
a X,fe, d rk2 vj )(d n Vi) + (d ri affl iM ) v { (d rki d Tk2 Vj ) 

_ a i,j;ki,k 2 

coefficient functions otf^ : 

[(d rki d rk2 d rk: vi){d ri Vj) + 3(d rki d rk2 vi)(d rk3 d n v j ) 

+ 3(d rki v l )(d rk2 d rk3 d ri v j )) (A4) 

At this point we see that for given n > l,m = there are n — 1 derivatives in the second 
term. The same is true for all terms except for the second in the contribution for n—1. But 
in general, i.e., except for n—1, these derivatives are distributed differently among the two 
v components in each term. Thus, ctf 1 ^ kn = const, follows for n > 2. Next, we will show 
that even ocf 1 ^ kn = for n > 2. This is immediately clear for odd n > 3 since respectively 
one particular term in the second line of (|A4|) has a unique derivational structure. For even 
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n > 4 there are respectively two terms, but with distinct derivational structure, leading 
again to cxfj?^ kn = 0. For n = 2 this only follows after evaluating the more complicated 
constraint 



The situation is different for n — 1, where the first and second term of ( |A2j ) lead to the 
constraint 



A lengthy evaluation yields the representation otf'^ = SijCt^ with some new functions 



fei • 



Now we turn to (|A1|) , the contribution for n = 0. The second term has to vanish, because 



it is the only one without any derivatives of v, which is only possible for off^ = const. The 



first term is of the same structure as the remaining second term of and the rhs. of (40), 
giving finally the constraint 

whose general solution 

atj = —ocf^ Tj — (z — l)ri + 7; with 7; independent of r (A5) 

we obtain after a few more calculational steps. 

coefficient functions otfl ■ 

aff\d L v j ){d n v i ) + (d n aff ] )vi(d L Vj) - aff ] (d L vi)(d ri Vj) (A6) 
coefficient functions afj 1 ^ : 

aff\dl Vj ){d nVi ) + (drATMd 2 ^) - a^Kdlv^Vj) + 2{d L v l ){d L d n v j )\ (A7) 

Similar to the contributions for n > 2,m = we find oif^ = for m > 2. For m = 1, on 
the other hand, ocf^ = const, holds because of the second term in fl^). And because of 
the other two terms we have even 

ctff 1 = o~ij/3i with Pl independent of r. (A8) 
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coefficient functions oq^ : 

( 9 r kl d LVj ) (d ri Vi) + (d ri agjjvz (d rki d L Vj ) 

- ^liiP^dLvdid^Vj) + (<9 rfci ^)(<9zA,^) + (d L vi)(d rki d ri Vj)] ( A9 ) 

This implies afj!^ kn = for n, m > 1. 
coefficient functions (3^°^ : 

i&fr^i)^) + (dr^lMd^Vj) - ^(d^drM) (A10) 

coefficient functions ^fj ■ 

&iM d Sk 2 v i){ d nVi) + (dn^l M )vi(d Ski d 8k Vj) 

~ tiiZkM^M 9 ^) + l(ds k Vi){d 8k2 d nVj )\ (All) 

As before we infer (3\ n f Q k \ kn = for n > 2. But for n = 1 /^'.^ = const, holds because of 
the second term in ( |A11|) , and furthermore we have 



Ajffci = $ijPki with independent of r. (A12) 
coefficient functions ftf*^ : 



- &(d Skl d ^( d ^ + (d Skl Vi)(d L d riVj ) + {dmW^drM)) (A13) 

It follows fij^fa kn = for n, m > 1. 

Combining the partial solutions (|A5|) -( |AT2D we are now able to write down the general 
solution (l2|) of (00). 
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FIGURES 

FIG. 1. Plane shear flow between two plates with distance L moving with relative velocity Ul 
FIG. 2. Illustration of definition (|3|) of the velocity increment w(r, i|xo, to) 
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